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We present new results on the cosmic star formation history in the SXDS-UDS field out to 
z = 1.6. We compile narrow -band data from the Subaru Telescope and the Visible and In- 
frared Survey Telescope for Astronomy (VISTA) in conjunction with broad-band data from 
the SXDS and UDS, to make a selection of 5725 emission-line galaxies in 12 redshift slices, 
■ spanning 10 Gyr of cosmic time. We determine photometric redshifts for the sample using 

£f~) ' 11-band photometry, and use a spectroscopically confirmed subset to fine tune the resultant 

redshift distribution. We use the maximum-likelihood technique to determine luminosity func- 
tions in each redshift slice and model the selection effects inherent in any narrow-band selec- 
| tion statistically, to obviate the retrospective corrections ordinarily required. The deep narrow 

. band data are sensitive to very low star formation rates (SFRs), and allow an accurate evalua- 

tion of the faint end slope of the Schechter function, a. We find that a is particularly sensitive 
to the assumed faintest broadband magnitude of a galaxy capable of hosting an emission line, 
and propose that this limit should be empirically motivated. For this analysis we base our 
K> ' threshold on the limiting observed equivalent widths of emission lines in the local Universe. 

5_j , We compute the characteristic SFR of galaxies in each redshift slice, and the integrated SFR 

density, psfr- We find our results to be in good agreement with the literature and parametrize 
the evolution of the SFR density as psfr oc (1 + z) 4 58 confirming a steep decline in star 
formation activity since z ~ 1.6. 

Key words: cosmology:observations - surveys - galaxies:evolution - galaxies:formation - 
galaxies:high-redshift - galaxies:luminosity functions. 



1 INTRODUCTION 

A key goal and surmountable challenge for observational cosmol- 
ogy today is to assemble a reliable picture of the formation and evo- 
lution of galaxies, and the build up of stellar mass in the Universe. 
The secular evolution of galaxies is a result of the formation of new 
stars, and so understanding the processes triggering and truncating 
star formation in addition to placing constraints on the evolution of 
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the star formation rate (SFR), in turn allows us to constrain theo- 
retical models of galaxy formation and evolution. A multitude of 
star formation indicators have been utilized in evaluating the SFR, 
however each of these is subject to its own biases and caveats, in 
addition to the limitations technology places on the redshift range 
where each method is applicable. 

Ultraviolet (UV) flux provides the most direct measurement 



of SFR by tracing hot young massive stars (e.g. Bunker et alj 
2004 lOuchi et alj |2004 rArnouts et al.ll2005L iBaldrv et all 12001 
Bouwens et al J I2010T) however UV measurements are subject to 
large dust corrections and so carry a large uncertainty in addition to 
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the danger of missing heavily obscured regions of star formation. 
The dust enshrouding young stars can itself be used to trace star 
forma tion, as it re-radiates the incident UV radiation in the infrared 
(e.g. IPerez-Gonzalez et aT1l2005t llxFloc'h et alj|2005l) . However, 
this method in turn struggles to detect some unobscured star forma- 
tion. Less commonly used methods include using radio emission to 
trace supernova remnants as sociated with short-lived massive stars 
(e.g. lMauch & Sadleill2007t) and tracing X-ray emission produced 
in X-ray binaries (of which one component is a massive O or B 
star). Alternatively, emission lines generated in H II regions can be 
used to trace either the recombination of atoms or the forbidden 
lines that arise in this low-density environment. Emission lines in 
the optical and UV are of course subject to dust corrections much 
as many other tracers, however Ha luminosity is particularly well 
calibrated for extinction effects. 

The evolution of the star formation rate density, psfr, is well 
studied, and a strong upward evolution is found with redshift out to 
z ~ 1 indicating that the p eak of star formation activi ty occurs at 
z > 1. The compilation o f Hop kins & Beaconj d2006h provides a 
broad overview of this relationship out to z ~ 6 but demonstrates 
that, even within the relatively well-constrained epoch, estimates of 
SFR evolution still vary according to indicator used. For instance 
the evolution of psfr is often parametrized as psfr. oc (1 + z)"' , 
and in one of th e earliest studies (using a rest frame UV colour 
to measure SFR) iLillv et alj d 19961) found 7 = 4. However, their 
sample was hampered b y a small survey are a , and the larger UV- 
sele cted samples of e.g. Cowie et"al] dl999l) . IWilson et alj d2002h 



and IPrescott et alj J2009 ) found much shallower slopes of 7 = 
1.5, 1.7 ± 1.0 a nd 2.5 ± 0.3 respective ly. Meanwhile, the FIR lumi- 
nosity found in lLe Floc'h et alj (12005) translates to a 7 = 3.2^2- 
Studies making use of continuum emission as a measure of the 
star formation activity of galaxies, whether it be UV emission, or 
its re-radiation in the infrared, allow us to quantify star formation in 
bright galaxies even down to low SFRs. The samples however are 
generally small, as sources not only require spectroscopy but must 
also be apparent in broad band filters in order to be detected — the 
imposed broad band flux limit means that samples in the optical 
become effectively stellar-mass limited, and therefore do not detect 
the lower mass population of galaxies responsible for an important 
contribution to the overall psfr. 

The use of narrow-band filters to detect emission-line galaxies 
provides a complementary technique to broad-band- selected sam- 
ples. A narrow-band survey selects objects based on the strength of 
line emission, and so while the equivalent width of emission lines 
comes into play, samples are close to SFR-limited, allowing the 
analysis of SFR over a much broader range of stellar masses. Re- 
cently, many authors have utilised this technique to compile large 
samples of star-forming galaxies and examine the evolution of the 
SFR as traced by eith er Ha luminosity (which scales directly with 
SFR; lKennicutdl 19921) or alternatively the [O III] and [O II] forbid- 
den lines which give a more loosely calibrated indication of the rate 
of star formation. 

Ha su rveys have been hug ely successful, from the pioneer- 
ing work of iBunker et alj d 19951) detecting only a handful of emit- 
ters, to the recent advances from the High Redshift Emis sion Line 
Survey (HiZELS: Icieach et alJ|20o3. ISobral et alj|2009l) that de- 
tected thousands of emission-line sources. Many authors have taken 
this approach, and the Ha lumin osity function is well-studied 
in a v ariety of environments (e.g. [ Gallego et al. 1995L Yan et alj 
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the Ha line beyond z = 0.4, strong emission lines blueward of Ha 
are required to study the higher redshift population using existing 
optical data s ets. The [O il] doublet for example has been traced 
to z = 1.6 bvlHogg et alj ill 9981). Ifflcks et alJ <2002h.lTeplitz et al 



2003) jAiild et dj {2003hjDrozdovskv et alj|2005UGilbank et al 



2010) andlLv et al] fe012h 



iFuiita et alj 12003k iKodama et alj 12004 



Umeda et alj 12004 IShiova et alj 120081 and IVillar et alj 120081 to 
name but a few). Since infrared spectroscopy is required to trace 



iLv et alj d2012h incorporated the [O III] tracer as well and 
derived luminosity functions in 11 redshift slices in the Subaru 
Deep Field based on Ha, [O III] and [O II] luminosities. Fitting 
a Schechter function to their data, they found strong evolution in 
the faint-end slope, a, and the characteristic number density <f>*, 
but little evolution in the characteristic luminosity of galaxies, L*. 
They interpreted these results as a comparatively stronger evolution 
of the low-luminosity galaxy population than those with higher in- 
trinsic luminosities. They found a steep evolution in psfr, in good 
ag reement with previ ous results, a trend confirmed by the work 
of ISobral et al.l d20 13h who compiled the first view of star forma- 
tion frorn_j_^0.4out to z — 2.2 using only Ha-selected galax- 
ies. ISobral et alJd2013l) parametrized the evolution in psfr as log 
Psfr = —2.1/(z + 1) suggesting a 7 ~ 3 across the redshift 
range < z < 1.6. Conversely however, they found no evolution 
of a, and set the parameter to —1.6 ± 0.08 for the last 11 Gyr, 
concluding that previous claims of evolution stemmed from het- 
erogeneous samples. This value of a was in good agreement with 
lHaves et alj d2010h who found a slope of a = -1.72 ± 0.20 at 
z — 2.2 for an Ha-selected sample. Even consistently selected 
samples however (i.e. those using the same emission line to derive 
each luminosity function) still show discrepancies in the value of 
a and its evoluti on with redshift (or lack of). The Ha emitters at 
z = 2.2 found in iTadakietalJ 4201 lh for example res ult in a slope 
of a = —1.37, much shallo wer than the studies of lHaves et alj 
(2010) and Sobral et al1 d2013l) and suggesting littl e evolution from 
high redshift to the present day. lHaves et al.ld2010h however, while 
in agreement with a steep value of a argued that their results did 
indeed confirm an evolution in a comp ared to the sha l lower slope 
of a = -1.35 at z = determined by iGallego etaLl dl995l) . Ev- 
idently, values of a across redshift remain unclear, and tentative 
claims of the parameter's evolution divide opinion further still. It 
is important therefore to better understand the factors which most 
affect the derivation of a, in order to place more meaningful con- 
straints on the faint end slope. 

This paper demonstrates how the maximum likelihood method 
can be used to determine luminosity functions, and the benefits 
of this approach over existing methods. Narrow-band selection is 
inherently sensitive to only a certain region of parameter space, 
and retrospective corrections are ordinarily applied in order to ac- 
count for the undetected objects. Here, we use a maximum like- 
lihood analysis to model the selection effects statistically and ro- 
bustly under assumptions about the underlying galaxy population. 
We explain in detail how we implement th is method, and co nfirm 
the technique by r eplicating the results of Sobr al et"al] d2012h and 
ISobraletal.ld2013h . 

Our luminosity functions are derived from narrow-band- 
selected emission-line galaxies in the Subaru/XMM-Newton Deep 
Field (SXDF), coincident with the Ultra Deep Survey (UDS) field 
of the UK Infrared Telescope (UKIRT)'s Infrared Deep Sky Sur- 
vey (UKIDSS). We use the Ha, [O III] and [O II] emission lines to 
determine field luminosity functions in 12 redshift slices, spanning 
10 Gyr of cosmic time, out to z ~ 1.6. We place constraints on 
the characteristic SFR in each redshift slice, and make estimates 
of Psfr for comparison to values in the literature. The outline 
of this paper is as follows. In Section [2] we summarize the data 



© 2012 RAS, MNRAS 000.[JJfl6] 



Evolution of Star Formation in the UDS out to z = 1.6 3 



used in this analysis, in Section [3] we discuss the source extrac- 
tion and selection of narrow-band emitters, and present the colour- 
magnitude plots of narrow-band selected objects. Section [4] dis- 
cusses our spectroscopically-confirmed subset, the star-galaxy sep- 
aration technique applied in order to remove late-type stellar con- 
taminants, our photometric redshift analysis and the resultant pho- 
tometric redshift distributions, and finally our method of redshift 
slice assignment for emitters selected in each filter. Section[5]details 
our application of the maximum likelihood technique for determin- 
ing luminosity functions, and discusses in detail the factors affect- 
ing the narrow-band selection technique. We describe also how this 
analysis carefully models the selection effects ordinarily corrected 
for retrospectively. In Section [6] we present luminosity functions 
(LFs) and derive the characteristic SFR and psfr in each redshift 
slice for comparison to literature results. Section [TJdiscusses the 
factors affecting our results and the effect of varying assumptions 
about the underlying galaxy population on the resultant values of 
a. We summarize our conclusions in Section|8] 

An Ho = 70 kms -1 Mpc -1 , SI m =0.3 and Qa =0.7 cosmology 
is assumed throughout, and all magnitudes are in the AB system. 
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Figure 1. Wavelength coverage and transmission functions for all selection 
filters. Four narrow-band filters on Subaru, plus two on VISTA (bold black 
lines) together with broad-band filters at corresponding wavelengths also 
from the Subaru Telescope (dashed lines) are shown. 



2 DATA 

The data used (summarized in Table[T) are taken from the SXDF- 
UDS field in order to make use of the deep optical data from the 
Subaru/XMM-Newton Deep Survey (SXDS), in addition to deep 
infrared imaging from UKI DSS UDS. 

The SXDS jFurusawa et ai]|2008l) covers 1.3 square degrees 
and delivers deep optical imaging in the B, V, R, i and z bands in 
five overlapping pointings of Suprime-Cam. 

UKIDSS, c omprised of 5 sub sur veys of various areal cover- 
ages and depths JLawrence e t al. 2007) targeted the SXD F with the 
deepest of these surveys, the UDS dFoucaud et al .120071 Almaini et 
al., [in prep]) covering 0.77 sq degrees. Here we utilise the Eighth 
Data Release (DR8) from this survey (details in TableQJ. Addition- 
ally we incorporate u band data across the entire field from Mega- 
cam on the Canada France Hawaii Telescope (CFHT), and Spitzer 
IRAC imaging from the Spitzer-UDS (SpUDS) survey in channels 
1 and 2. 

Narrow-Band imaging was obtained from Suprime-Cam in 
four narrow-band filters, and smoothed using a Gaussian kernel to 
recover the same point spread function (PSF) as the broad-band op- 
tical data fr om Subaru (0.82 ar csec onds). Details of th ese data can 
be found in lOuchi et ail d2008l) and lOuchi et"al] j2009h . 

In addition, we make use of two narrow-band filters on the 
Visible and Infrared Survey Telescope for Astronomy (VISTA ). 
VISTAs visible-infrared camera (VIRCAM; balton et al.ll2006h . 
is a wide-field near-IR camera consisting of sixteen 2048 x 2048 
Raytheon VIRGO HgCdTe arrays distributed over the focal plane 
of VISTA. The paw-print of the 16 detectors provides a non- 
contiguous instantaneous field of view of 0.6 deg 2 . Half of detector 
16 is affected by a time- varying quantum efficiency which makes 
flat fi elding that detector extremely challenging (see Ijarvis et al] 
|2013| for further details). The NB980 and NB990 filters (Orr, Wal- 
lace and Dalton 2008) share a slot on VIRCAM — each covering 
eight detectors, and so to obtain a full pawprint from either of these 
filters the survey area has to be observed twice, with the instrument 
rotated 180 degrees. Consequently, the two VISTA NB images suf- 
fer from mismatched seeing in the top and bottom halves of the 
images, and so we smooth the data to the PSF of the worst-seeing 
parts of the VISTA pawprint (1.42 arcseconds). 



In order to select a reliable sample with homogeneously de- 
termined photometric redshifts, we use only the region of sky cov- 
ered by both SXDS and UDS imaging so as to provide 1 1 bands of 
photometry for every detection. The usable overlap of these deep 
broad-band surveys and the Subaru and VISTA narrow-band imag- 
ing data amounts to 0.63 sq degrees and 0.38 sq degrees respec- 
tively. The 6 narrow-band filters used for selection are shown to- 
gether with Suprime-Cam broad-band filters across the correspond- 
ing wavelength ranges in Figure [T] 



3 SAMPLE CONSTRUCTION 

We use the Source extraction software SEXTRACTOR 
l lBertin & Arnoutslll996l) to detect objects in each of the narrow- 
band (NB) images. SEXTRACTOR works by constructing a 
background map and then looking for groups of connected pixels 
above a threshold calculated from the root mean square (RMS) 
noise in the image. We require a group of 10 connected pixels 
at a significance of 1.5 times the background RMS to confirm 
a detection. We then estimate the narrow-band image depths by 
placing 10 000 randomly positioned 2-arcsecond apertures on 
each image and iteratively fitting a Gaussian to the sigma-clipped 
histogram of resultant fluxes. The catalogue is then cut to demand 
an NB-detection above 5a relative to our estimate of the sky noise 
in order for a source to be selected as a possible line emitter. 

We make a very conservative area cut, dicarding first the offi- 
cial SXDS masked areas, and then visually inspecting the narrow- 
band images (with object detections overlaid) to discard additional 
areas of noisy coverage towards the edges of the CCD and then ar- 
tifacts, areas around bright stars and CCD bleeds in the broad-band 
images. The areas removed from the study as a result of the visual 
inspection varied slightly from filter to filter. 

3.1 Subaru Selection 

For objects selected in one of the four Subaru NB filters, the sample 
was culled to the 5a limiting magnitude of the shallowest SXDS 
field. We next consider the interpolated broad-band colour at the 
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Table 1. Summary of photometric data in all narrow and broad-band filters from Subaru, VISTA, CFHT, UKIRT and Spitzer. 



Telescope; Instrument 


Filter 


Filter Type 


Central wavelength 


r W rllvl 


5 sigma limit (2" aperture) 


Citation 


Subaru: Suprime-Cam 


NB503 


NB 


5029A 


/J A 


24.91 


oXD^, Uucni et al. (zUUo) 


Subaru: Suprime-Cam 


NB570 


NB 


5703A 


68A 


24.52 


oXDc>, Uucni et al. (zUUo) 


Subaru: Suprime-Cam 


NB816 


NB 


8150A 


1 19A 


25.65 


SXDS; Uucni et al. (20Uo) 


Subaru: Suprime-Cam 


NB921 


NB 


V to JA 


1 JlA 


25.38 


aXDo, Uucni et al. (zUUv) 


Viol A. V1KL.A1V1 




MB 
IN D 


y/otA 


ins A 


Tl 1-7 

Zj. t / 


Jarvis et al., (in prep) 


VTSTA- VTRCAM 


NB990 


NB 


9909A 


105A 


23.54 


Tnrvn; pt nl fin nrpn \ 


CFHT: Megacam 


u 


BB 


3740A 


500A 


26.72 


Foucaud et al., (in prep) 


Subaru: Suprime-Cam 


B 


BB 


4473A 


1079A 


27.39 


Furasawa et al. (2008) 


Subaru: Suprime-Cam 


V 


BB 


5482A 


984A 


27.10 


Furusawa et al. (2008) 


Subaru: Suprime-Cam 


R 


BB 


6531A 


1160A 


26.85 


Furasawa et al. (2008) 


Subaru: Suprime-Cam 


i 


BB 


7695A 


1543A 


26.66 


Furasawa et al. (2008) 


Subaru: Suprime-Cam 


z 


BB 


9 149 A 


1384A 


25.95 


Furasawa et al. (2008) 


UKIRT: WFCAM 


J 


BB 


12500A 


1570A 


24.98 


UKIDSS UDS; Almaini et al., (in prep) 


UKIRT: WFCAM 


H 


BB 


16500A 


2910A 


24.27 


UKIDSS UDS; Almaini et al., (in prep) 


UKIRT: WFCAM 


K 


BB 


22000A 


3530A 


24.59 


UKIDSS UDS; Almaini et al., (in prep) 


Spitzer: IRAC 


chl 


BB 


35500A 


741 lA 


23.55 


Spitzer SpUDS: PI Dunlop 


Spitzer: IRAC 


ch2 


BB 


44900A 


10072A 


22.88 


Spitzer SpUDS: PI Dunlop 



wavelength of the selection filter as a measure of the continuum 
level for each object. We calculate this continuum (IntNB con t) by 
interpolating between the two broad-bands adjacent to the NB as- 
suming that continuum flux follows a power law: 



IntNB cont = BB1 - (BB1 - BB2) ^(WW , 

log(ABB2/ABBl) 



(1) 



As filter NB921 lies centrally in the z-band, we calculate 
colour in this filter relative to the z-band only. Broad-band magni- 
tudes are obtained by running SEXTRACTOR in dual image mode, 
and colours are measured in 2 arcsecond apertures. Objects are ex- 
pected to scatter around zero in (IntNB cont — NB), with emission- 
line sources showing a positive colour. In practice, the median 
colour of NB-selected sources is slightly offset from zero, due to 
the mismatch of interpolated broad-band magnitude and the true 
continuum magnitude at the wavelength of the NB. Furthermore, as 
we select objects across a broad range of NB magnitudes, we are es- 
sentially sampling different redshifts (in addition to different parts 
of the galaxy population) and so with changing NB-magnitude, we 
introduce a systematic error that gives rise to a gradient in the me- 
dian colour of objects (typically of order ~ 0.015 mag). We remove 
the colour gradient by calculating the sigma-clipped median colour 
of sources as a linear function of NB mag, and subtracting this so 
that the resultant locus of colours lies along zero. 

Potential emission-line objects must display an NB-excess of 
at least 3a relative to the sigma clipped median colour and 3a rel- 
ative to the scatter from photometric uncertainty (in the relevant 
SXDS field) to be selected as an emitter. Photometric scatter is de- 
termined from the image depth estimates described above. Figure 
[2]shows colour-magnitude diagrams for each of the NB selections, 
with green lines showing the lines of 1, 2 and 3 a colour excess, and 
the red lines depicting the 1, 2 and 3 a photometric error trumpets. 

We assume our observations are background limited, as even 
at the bright magnitudes where photon shot noise comes into play, 
colour scatter in (IntNB con t-NB) is determined empirically and so 
the requirement of a 3a excess takes this scatter into account. 



3.2 VISTA Selection 

For selection in one of the two VISTA filters (seen alongside the 
Subaru filters in Figure [T}, detections are made and NB-excesses 
evaluated following broadly the same method developed with the 
Subaru data. The only variation in the selection process is the man- 
ner in which we evaluate the continuum contribution to the NB 
measurement. The two VISTA NBs are close in wavelength, and so 
they sample objects at (effectively) the same point in their emission 
spectra; measured fluxes should therefore only differ significantly 
when an emission line is present in one of the two filters. This al- 
lows us to detect objects in one NB filter, and use photometry from 
the other as a measure of continuum flux for each object. 

The double-NB selection technique circumvents the error in- 
troduced through inaccurate determination of interpolated broad- 
band colour (the main source of scatter introduced into the NB- 
broad-band colour) and so we anticipate the locus of VISTA- 
selected objects to be much tighter, allowing for a much less conser- 
vative 3a cut in colour scatter, i.e. a smaller colour-excess is now 
robust at the 3a level. Colour-magnitude diagrams for these two 
NBs are also shown in Figure[2]demonstrating the anticipated nar- 
row locus, however the far shallower depth of the VISTA imaging 
results in a shallower sampling of the star-forming population. 



4 REDSHIFTS 

Our narrow-band selection, spanning six NB filters from 5029 A to 
9909 A is sensitive to all the major diagnostic emission lines in the 
region from [O II] to Ha + [N II], and samples redshift slices from 
z=0.14 ([O III] in the NB570 filter) to z=1.6 where [O II] falls in to 
the VISTA filters. With more than 7000 objects meeting the selec- 
tion criteria, it is essential that we are able to reliably discriminate 
between the different lines falling in each filter while employing 
a technique suitable for dealing with such a volume of data. The 
combination of the deep SXDS-UDS imaging, and complementary 
u band and IRAC data, provides 11 broad bands of photometry 
available for every detection, allowing reliable determinations of 
objects' spectral energy distributions (SEDs) and consequently ac- 
curate photometric redshifts. We have obtained spectra for a small 
subset of objects (discussed in Section l4~3t and use the redshifts 
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Figure 2. The colour-magnitude diagrams for the four Subaru narrow-band filters: NB503, NB570, NB816 and NB921, and the two VISTA narrow-band 
filters: NB980 and NB990. Detections are restricted to the 5<r limiting NB-magnitude in the shallowest SXDS field. The solid green line depicts the median 
colour of sources (corrected for gradient and offset from zero — see text for details), dashed green lines show the 1, 2, and 3a limits on intrinsic colour scatter. 
Solid red lines depict the 1, 2, and 3<r photometric errors respectively for the shallowest region of the five SXDS fields. Solid blue circles represent the entire 
sample of NB-excess-selected objects. 
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Figure 3. BzK colour-colour diagrams for the two reddest filters in our Subaru NB-Excess selection; NB816 shown on the left and NB921 on the right. All 
objects falling below the solid black line are classed as stars according to Daddi et al. ( 2004). The subset of these objects that we identify as stars (and therefore 
eliminate from our analysis) are highlighted by red stars, yellow diamonds represent confirmed stars in the UDS. See main text for details. 



derived as a training set for the photometric redshift code "EAzY" 
(Brammer, Van Dokkum and Coppi 2008). 



4.1 Photometry 

Photometry for the full selection of NB-excess selected objects is 
computed in all 1 1 broad bands using the IRAF task phot in 2 arc- 
second apertures. Astrometric corrections between Suprime-Cam 
data and the other im ages were performed using the method of 
Simps on" et alj J2012h . All Suprime-Cam imaging was reduced to 
provide a uniform PSF of 0.82 arcsec, and the J, H and K band 
images from WFCAM are consistent with this to within 0.05 arc- 
sec, and so no corrections are applied to the measured photome- 
try. B and K band images were smoothed using a Gaussian ker- 
nel to recover the same PSF as the u and IRAC images respec- 
tively, and fluxes from the u band image and the two IRAC images 
were then scaled according to the ratio of fluxes in the unsmoothed 
and smoothed images. All broad-band photometry i s corrected for 
Galac tic extinc tion of A v = 0.070 magnitudes dSchlegel et alJ 
1 1998b using the lCardelli. Clayton, & Mathisl dl989h reddening law. 



4.2 Star-Galaxy Separation 

Our narrow-band-excess technique is inherently susceptible to the 
spurious selection of any object with a spectrum not well described 
by a power law. This leads to contamination of the sample by late- 
type stars which are known to exhibit very red, complex SEDs due 
to the many absorption features present in their spec tra. We employ 
the two-colour BzK technique o fDaddi etal.U2004h to identify and 
eliminate these stars from the red filters susceptible to this contami- 
nation. For the Subaru NB filters this applies to NB816 and NB921. 
Daddi et al., (2004) derive Equation[2] and class all objects below 
this line as stars. 



(z - K) < 0.3 (B - z) - 0.5 (2) 
As we are contaminated only by late-type, red, stars, we em- 



Table 2. Detections per NB filter. Excess objects are defined as objects pre- 
senting an NB-excess above 3<r in intrinsic colour scatter and photometric 
error. Galaxies are those objects that remain following the elimination of 
stars via our star-galaxy separation technique (see Section l4~2l . 

Filter Detections Excess Objects Galaxies 



NB503 


43781 


4368 


4368 


NB570 


37846 


1749 


1749 


NB816 


97382 


3816 


3595 


NB921 


109126 


5395 


5360 


NB980 


16193 


119 


119 


NB990 


17755 


191 


191 



ploy the additional criterion that objects must have B-z >2.5 to be 
classed as stars and eliminated from the selection. 

The BzK colours of our sample are s hown in Figure [3] The 
population of blue objects falling below the Daddi et al. (2004J) BzK 
solid black line arises from scatter as a consequence of very faint K 
magnitudes (> 26). We eliminate 256 stars in total, amounting to 2 
per cent of the NB-excess sample. 

Despite the red wavelength of the two VISTA filters, only 
NB990 is in danger of the spurious selection of stars due to the 
point in a late type star spectrum it samples. Running the NB990 
selection through our BzK criteria however eliminates nothing from 
the small sample. Table [2] summarizes the number of detections in 
each filter plus numbers remaining in the sample after quantifying 
NB-excess and removing stars from the sample. 



4.3 Spectroscopic Subset Analysis 

With a preliminary selection of NB-excess objects in the central 
field of the SXDS, we compiled a target list for follow-up spec- 
troscopy. Observations were carried out with the Inamori Magel- 
lan Areal Camera Spectrograph (IMACS) instrument on the Walter 
Baade Telescope on the dates of 201 1 Jan 03-04, using 0.8 arcsec- 
ond slits and the 200 lines per mm grism (covering the wavelength 
range AA = 0.39 - 1.00/xm, giving AA ~ 8 A pix -1 at FWHM). 
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Figure 4. Comparison of spectroscopic and photometric redshifts for the 
subset of objects with spectroscopic follow-up, and spectroscopic redshifts 
of the highest quality flag. The bottom panel shows £ p hot 

against z spec for 

the optimized set of EAzY parameters. The top panel shows the distribution 
of Az around zero. Cyan points represent our Magellan spectroscopic data, 
and black points represent all other available spectra in the UDS. 



The total exposure time was 335 minutes, and seeing varied from 
0.6 — 0.9 arcseconds. 

The spectra were reduced using the COSMOS reduction 
pipeline befor e extracting redshift in formation using the runz red- 
shifting code dSaunders et al.l 12004) . We determine reliable red- 
shifts for 129 of the sample, and assign 108 of these the highest 
quality flag. The remainder of the subset is made up of low signal 
to noise detections (where emission-line flux does not in fact meet 
our final criteria), spectra displaying a single unidentifiable line, 
spectra affected by cross-talk from other slits, and a small number 
of stars targeted to confirm our star-galaxy separation technique. 

The modest but reliable subset of spectroscopic redshifts con- 
firms that we are indeed selecting emission-line galaxies at the red- 
shifts anticipated, and serves as a reference to determine the suc- 
cess of our photometric redshift technique applied to the rest of the 
sample. 



4.4 Photometric Redshift Analysis 

The default parameters of EAzY provide photometric redshifts in 
good agreement with the distribution expected for emission-line 
objects selected in these filters; however, to optimize the SED- 
fitting, we modify this default parameter set so as to minimize 
onmad for the spectroscopic subset (Equation [3] B rammer et al.l 
2008), where 



cnmad 



1. 



x median 



Az — median(A2 



1 + z B 



(3) 



The number of catastrophic outliers (> 5ctnmad) for definite 
spectroscopic redshifts (of the highest quality flag) equates to 5 per 
cent of the spectroscopic subset, with a value of ctnmad = 0.026. 
Including objects with less certain spectroscopic redshifts raises 
this to 9 per cent. While photometric redshifts are commonplace for 
large survey data sets, particularly for A"-selected (i.e. stellar-mass- 
limited) samples, photometry is frequently restricted to broad- 



bands only. Here we possess additional information that all objects 
remaining after the selection process must display a narrow-band 
excess at the wavelength of the filter. 

While naively, we might interpret the presence of a narrow- 
band excess to mean that all sources must possess a major emission 
line redshifted to the wavelength of the selection filter, in prac- 
tice the problem proves to be more complex. We design a pure 
emission-line template to investigate the effect of including narrow- 
band data in addition to broad-band photometry, before fine tuning 
the setup parameters. Our emission line template consists of the 
major emission lines Lya, [O II], H/3, [O III], Ha and [N II] with 
negligible continuum flux. However, including this template now 
means that the complexities of the emission-line ratios come into 
play, along with the nature of the narrow-band excess. While the 
objects do indeed lie in distinct redshift slices, narrow-band excess 
can arise from a vast number of emission lines. We find that the in- 
clusion of narrow-band data, along with an emission-line template, 
biases redshift determination towards the redshift slices associated 
with the emission lines in our template at the expense of a good fit 
to the broad-band photometry, thus adding redshift interlopers to 
our selection. The 1 1 bands of photometry are deep, and as EAzY 
is well equipped to determine redshifts based on this information 
alone we discard the narrow-band plus emission-line template setup 
at this stage. 

We find the best agreement between z p h ot and z Bpcc for the 
default template set, using broad-band data only (see Figure |4j. 
This way redshifts are assigned independently of their detection 
method, allowing the correct identification of spectral breaks and 
additional better fitting emission lines than those included in the 
emission-line template. For the comparison in Figure |4]we also in- 
corporate matches to the full set of spectra availab le in the UDS 
field. The larger SXDS-UDS spectroscopic sample (Simps on et alj 
I2012L and papers cited therein) includes many spectra that are the 
result of targeted follow-up of X-ray or radio sources, and is there- 
fore heavily biased towards Active Galactic Nuclei (AGN). Objects 
with spectra from the larger sample (depicted by black dots in Fig- 
ure [4) often appear at redshifts other than those typical for our se- 
lection, as they sample many objects with emission lines usually 
weak in star-forming galaxies. 



4.5 AGN Contamination 

We identify 8 (of the 129) ob j ects a s AGN according to the line- 
ratio diagrams of lLam areill el d2010b or based on the presence of 
[Ne III] emission. As we can use only our own spectroscopic 
follow-up from the Magellan telescopes to quantify AGN contami- 
nation in an unbiased manner, we infer a contamination fraction of 
6-39 percent. This correction is small and highly uncertain, and 
because AGN activity is undoubtedly associated with some star 
formation, removing these sources then leaves us susceptible to a 
possible over-correction and consequently an underestimate of star 
formation activity. We therefore choose not to correct our results for 
AGN contamination. 



4.6 Resultant Photometric Redshift Distribution 

EAzY provides an output probability density distribution P(z) for 
each object it evaluates. We use z m 2, the redshift marginalized over 
the posterior redshift probability distribution (also taking into ac- 
count the A"-band apparent magnitude), to determine the redshift 
distribution of our objects and discriminate between the emission 
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lines giving rise to narrow-band-excess. The resultant photometric 
redshift distribution can be seen in Figure [5] The advantage of the 
posterior redshift probability is two-fold: we combine the best esti- 
mate of redshift (using all SED templates) with a priori constraints 
on the likelihood of an object being at a particular redshift (given 
our knowledge of its colour and K magnitude prior), while taking 
into account the probability that the object lies at an alternative, 
less well-fit redshift. z m 2 effectively provides a weighted mean, 
however the resultant distribution of redshifts and that for the best- 
fitting redshift (z a ) are in fact very similar (traced by the two his- 
tograms in Figure [5}, which demonstrates that we are not altering 
the best-fitting redshifts dramatically through marginalisation and 
use of the prior. We note that the small number of catastrophic out- 
liers discussed in Section l4l4l produce doubly-peaked P(z) distri- 
butions, with z a often selecting a redshift of relatively small proba- 
bility compared to the probability at the value of z m 2- Similarly, the 
small peaks of objects lying at z a =0 are the result of objects miss- 
ing a red band of photometry and producing a small peak at z a =0 
in their P(z) distributions. This lends further support to the advan- 
tages of z m 2 which in these cases avoids the small z a =0 spike via 
the weighted mean. Figuref5]shows that we are selecting emission- 
line galaxies at the anticipated redshifts and highlights the peaks in 
the distribution associated with Ha, [O III] and [O II], in addition 
to Lya in the two bluest filters (Ouchi et al. 2008]; not included in 
this analysis). 

The lowest panel of Figure [5] shows the combined data sets 
of NB980 and NB990. As these two filters are very close in wave- 
length (A A c ~120A) we choose to combine the two data sets to 
combat the small number of detections (with a 5<r detection cutoff 
in the shallowest filter; NB980). The NB980 and NB990 data sets 
are evaluated separately, and then combined before assigning ob- 
jects to a redshift slice. The combined data set (hereafter NB985) 
is treated as a single selection from here on. 

4.7 Redshift Slice Assignment 

Before assigning objects to a particular redshift slice we consider 
the stack of P(z) values for each NB-selection of objects (seen in 
Figure [6j. We bin objects into the redshift slices associated with 
the peaks in the data, and assign the upper and lower limits of the 
bins according to the minima in the stacked P(z) distribution. For 
the VISTA selection where the minima of the distribution are less 
apparent, the widths of the NB921 redshift windows are applied, 
centred on the peaks of the VISTA data. Table[3]summarizes the fi- 
nal numbers of objects in each redshift slice and the redshift ranges 
over which each peak was defined. 

4.8 Line fluxes 

Narrow-Band selections are by necessity made in an observable pa- 
rameter space, here based on line flux. Apparent aperture line fluxes 
(hereafter: line fluxes) for objects in each redshift slice are derived 
from the initial SEXTRACTOR detection measurements. The en- 
tirety of the NB-excess is assumed to arise from an emission line 
and so line flux (/iJy) is computed as Equation|4] and converted to 
Wm" 2 . 

Line flux = Huxnb — flux cont (4) 

where Auxnb is the flux in the narrow-band, and flux con t is the 
continuum contribution to that narrow-band flux. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 



Figure 5. Photometric redshift distributions for the 4 Subaru narrow -band 
niters and the combined distributions from VISTA's NB980 and NB990 fil- 
ters. The dark grey filled histogram represents z m 2, the redshift marginal- 
ized over the posterior redshift probability distribution, and the light grey 
hatched histogram traces z a , the best fit redshift. Panels top to bottom: 
NB503, NB570. NB816, NB921 and NB980+NB990 in redshift bins of 
width Az= 0.1. Vertical lines indicate the redshifts at which different emis- 
sion lines would be detected in the various filters: triple-dot-dashed lines 
show Ha (dark red), dashed lines the [O III] doublet and H/3 (green), and 
dot-dashed lines [O II] in blue. Filters NB503 and NB570 also show peaks 
at z~3 where we detect Lya (dotted purple line). 



Table 3. Number of objects assigned to each redshift slice. The redshift 
range quoted refers to the limits in z m 2 according to which objects were 
binned into redshift slices. These limits are taken from the minima in Figure 
[6] with the exception of the VISTA data for which the NB921 limits were 
transposed to the appropriate redshift. The number of objects per slice is the 
number of real objects assigned to a redshift slice from our catalogue, that 
were used to derive a luminosity function for this redshift range. 



Filter 


Redshift 


Line 


Objects 


NB503 


0.10 < 0.35 < 0.50 


[Oil] 


142 


NB570 


0.00 < 0.14 < 0.30 


torn] 


42 


NB570 


0.30 < 0.53 < 0.70 


[Oil] 


96 


NB816 


0.00 < 0.25 < 0.35 


Ha 


152 


NB816 


0.35 < 0.63 < 0.80 


[O III] 


985 


NB816 


0.80 < 1.19 < 1.50 


[Oil] 


1013 


NB921 


0.00 < 0.40 < 0.50 


Ha 


279 


NB921 


0.50 < 0.83 < 1.10 


[O III] 


930 


NB921 


1.10 < 1.46 < 1.90 


[Oil] 


2204 


NB985 


0.00 < 0.50 < 0.60 


Ha 


153 


NB985 


0.60 < 0.99 < 1.29 


[Oni] 


33 


NB985 


1.29 < 1.64 < 2.08 


[Oil] 


23 
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Ha+ [Nil] 





[Oil] 



0.36 0.38 0.40 0.42 0.44 



0.75 0.80 0.85 0.90 95 



1 .40 1 .42 1 .44 1 .46 1 .48 1 .50 1 .52 



Figure 7. Line transmission for NB921. Panels from left to right show the effective line transmission with redshift for Ho, [O III] and [O II] respectively, as a 
synthetic spectrum is propagated through the NB921 filter. 




0.0 0.5 1.0 Z hnt 1.5 2.0 2.5 



Figure 6. Stacked probability density distributions for each NB-selection 
of objects, shown over the redshift range relevant for this study. Red (Ha), 
green ([O III]) and blue ([O II]) shaded regions represent the redshift ranges 
assigned to a slice for a particular filter/line combination. 



5 MAXIMUM LIKELIHOOD LUMINOSITY FUNCTIONS 

We em ploy the maximum likelihood method of Marsh all et alj 
( 1983) in order to fit luminosity functions to the data. We assume 
the distribution of line luminosities in a given redshift slice will 
follow a Schechter function, given in log form by Equation [5] 

cj> (L) dlogL = InlO (j>* ( jp J e~ (L/L 5 dlogi (5) 

where <j>* , L* and a are the characteristic number density, charac- 
teristic luminosity, and the gradient of the faint-end slope respec- 
tively dSchechtenll976t) . We assume negligible cosmic evolution 
across the small redshift interval per slice. Next we consider how 



the distribution of true line luminosities maps onto our observed 
distribution of line fluxes. 

Splitting the line flux range where we are sensitive to lines into 
bins small enough to expect no more than 1 object per bin, we can 
write the likelihood of finding an object in bins Fj and no objects 
in bins Fj, as Equation|6]for a given luminosity function: 

A = n*(*D dlo gj Fe-* (F ' )dlogF e -*W dl °^ (6) 

where fy(Fi) is the probability of detecting an object with line flux 
between F and 10 dlogF F. This simplifies to Equation|7] where Ffc 
is the product over all bins: 

A = I] dl °S F II e~* (Ffc)dlogF (7) 

Fi F k 

We then define the likelihood function as S = — 21nA and 
minimize this (Equation [8} in order to determine the parameters of 
the model ((f)* , L*, a) that best reproduce the observed data. 

S= -2^1n*(J , ! ) + 2^*(F)dlogF (8) 

To obtain $ (F) (the observed line luminosity function) from 
</> (L) (the true line luminosity function) we must consider the se- 
lection effects imposed in observing the data and statistically model 
these perturbations to cf> (L) before minimising the likelihood func- 
tion. 

5.1 Treatment of the filter profile 

To produce a LF, we must first assess the NB filter response for the 
emission line of interest as a function of redshift. A given line flux 
may arise from a line near the centre of of a detection filter, or an 
intrinsically brighter line lying closer to the edge of the filter where 
transmission is lower. 

We define T(z), the relative transmission with redshift, to de- 
termine the redshift range to which each filter/line combination is 
sensitive, accounting for filter transmission in addition to contribu- 
tions from other lines. T(z) is calculated by propagating a synthetic 
spectrum through the filter transmission curve (an example of the 
resultant line transmission is seen in Figure [7), using standard line 



© 2012 RAS, MNRAS 000. [01161 



10 Alyssa B. Drake et al. 



ratios ([N II]/Hq=0.33, [O IH]/H/ 3 = 3, [O II] in the low density 
limit; lOsterbrock & Fe rland 2006) normalized to the line of inter- 
est. F, our line flux, now becomes: 

where L / Aivdj (z) is the true flux, convolved with filter trans- 
mission T(z). ^(F), the probability of seeing a line of flux F, now 
becomes: 

*(F)dl 0g F = / (^f) (f ) d,dlogL (10) 

where dV/dz is the co-moving volume interval in our field per 
unit redshift, and the integral is performed over the redshift range 
to which the filter is sensitive. 

5.2 Line detectability 

The detection of an emission line in our study requires it to reside 
in an observable part of parameter space in the NB magnitude - 
colour-excess plane, i.e. the region within the three dotted lines in 
Figure [8] The detection fraction as a function of line flux, is then 
given by Equation QT] 



selection criteria. For an object to be included in our analysis it 
must display a colour excess of at least 3a significance. At bright 
magnitudes (where intrinsic colour-scatter dominates uncertainty) 
detections are limited by an equivalent width limit (dotted hori- 
zontal line in Figure [8} . Therefore, for a given line flux, only host 
galaxies fainter than a particular broad-band magnitude will result 
in an equivalent width large enough to be selected in our sample. 

For very bright lines however (such as that traced on the left 
hand side of Figure [8} the limiting factor becomes the saturation 
limit for point sources in the broad-band imaging used to determine 
narrow-band excess. Very bright lines (of flux greater than this) 
will not be selected. The value of this limit varies from filter to 
filter, and we define each limit as the point at which the locus of 
stellar objects crosses the 3a color-excess line. Brighter than this 
limit, emission-line objects could be indistinguishable from stellar 
objects. The broad-band saturation limits for each NB filter are as 
follows: NB503=19.0, NB570=19.1, NB816=19.25, NB921=19.2, 
NB980/NB990=16.0. We note that the effect of this limit is small, 
indeed removing this correction from the analysis altogether and 
integrating Equation[8]to infinite flux has no effect above z ~ 0.6, 
and very little effect below this. An example of a particularly bright 
line (of 60 a significance) is shown in Figure[8]depicted by the line 
on the far left. This demonstrates how only lines of very large line 
flux are limited by this saturation limit (point A2 in Figure |8j as 
opposed to the EW limit (Al). 



/det(F) 



''faint lin 
'-bright iii 



/>BB(mBB)dmBB 



?>BB(?riBB)dmBB 



(11) 



where <^bb("1bb) is a Schechter fit to the interpolated broad- 
band magnitudes of all galaxies within a particular filter's redshift 
range, and is der i ved from the photometric redshift catalogue of 
Grutz bauch et alj d201 lh . The numerator of /a et denotes the sam- 
ple of galaxies our technique is capable of detecting (thick red 
lines in Figure [8}, and the denominator denotes the entire popu- 
lation of galaxies capable of hosting the line down to some in- 
tegration limit, mf a i n t. mfaint is the assumed faintest magnitude 
of a galaxy capable of hosting a particular emission line luminos- 
ity. Different authors set this limit differently and somewhat arbi- 
trarily. Here we adopt a limit based on the rest frame equivalent 
width of the line o f interest seen in the local Universe, and set 
EW mfaint = 100A jKennicutl|[l992T) . The resultant Schechter fit 
to our observed line fluxes is sensitive to the integration limit of the 
denominator (mf a i nt ). The effect of varying this limit is discussed 
in Section |T2l 

The limits of integration for the numerator are somewhat more 
complex. The detection of a galaxy of particular broad-band mag- 
nitude for a line of particular line flux is limited at the bright end by 
the fainter of 2 broad-band magnitudes: either that corresponding 
to the EW limit, or that corresponding to the bright NB limit. At the 
faint end, broad-band magnitude is limited by the brighter of two 
broad-band magnitudes: that corresponding to the 5a NB detection 
limit, or that corresponding to the faintest broad-band magnitude 
capable of hosting a line of this line flux (i.e. mf a i nt ). Figure[8]il- 
lustrates the location of the bright and faint broad-band limits in 
relation to lines of constant line flux. 



5.2.2 Broad-band magnitude faint limit 

The faint limit for the majority of lines is set by mf a i nt , i.e. the true 
limiting magnitude of the population of galaxies capable of hosting 
that line, as detailed above. 

For exceptionally faint lines however, the detection limit of 
our catalogue comes into play. For an object to be detected its 
broad-band magnitude must be bright enough that when added to 
its line flux its NB magnitude is above the 5a detection thresh- 
old, and so for very faint lines this constraint sets the faint limit on 
broad-band magnitude. Point B 1 represents this limit for a particu- 
lar line flux in Figure [8] 

Equation[TT]then becomes Equationll2l 

rmin(m 5(7 ,m faint ) , / n , 

m \ <PBB(m.BB)aro,BB 

r I n \ Jmax l m EWlim> m NBbright.) v ' , , ,,, 

J_co 0BB(m B B)dm B B 
and 9(F) now takes its final form: 



(13) 



*(F)dlogF = / det (F) J <P (j^fj ^dzdlogK 

We note that the pre sence of a strong correlat ion between 
SFR and stellar mass (e.g. lNoeskdl2008l however see lSobral et ail 
1 201 ll for an alternative view) would mean that low-luminosity lines 
would be less likely to be found in more massive galaxies. We will 
address this issue and its effect on line detectability in paper 2, 
where we will study luminosity functions as a function of stellar 



5.2.1 Broad-band magnitude bright limit 

For the majority of lines the bright limit is set by mEwiim (depicted 
as point Al in Figure [8} and can be understood by considering our 



5.3 Detection Completeness 

Finally, we assess the detection completeness of the narrow-band 
selection catalogue retrieved via SEXTRACTOR. The completeness 
of detections is dependent on the depth of the detection image and 
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— Line of constant line flux 

Observable line fluxes 

A1 - BB bright limit for most lines 
A2 - BB bright limit for very bright lines 
B1 - BB faint limit for very faint lines 




Narrow-band magnitude 

Figure 8. Lines of constant line flux in the NBmag - (colour-excess) plane. 
The dotted lines represent the limiting equivalent width for an emission 
line to meet the selection criteria, and the faintest broad-band magnitude 
that when combined with a given line flux will meet the 5 a narrow-band 
threshold. A and B represent points on a line of constant line flux that corre- 
spond to the bright and faint broad-band limiting magnitudes respectively. 
Thick red lines highlight the portion of each line of constant line flux that is 
detectable in our sample. 



6 RESULTS 

We derive luminosity functions to examine the evolution of the Ha, 
[O III] and [O II] luminosity functions. The resultant Schechter fits 
provide LFs in 12 redshift slices between 0.14 and 1.64. Using the 
parametrizations of these LFs we can derive SFRs and compute 
Psfr to examine the evolution of star formation with redshift. 



6.1 Luminosity Functions 

The derived luminosity functions are presented in Tableland Fig- 
ure [9] Volumes quoted are representative of the volume covered 
where filter transmission is greater than 0.5. Observed (aperture- 
corrected) values are presented, followed by the extinction- 
corr ected value of L * , assuming 1 magnitu de of extinction a t 
Ha dKennicuttll 19921 ; 1998) and as appl ied i nlFuiita et ail |2003|) . 
IPascual et al.1 ( |2005|) . ISobral et all d2012h and lSobral etail llOll) . 
Figure [9] presents the best fitting LFs for the data at each redshift 
with the maximum likelihood fit shown in bold coloured lines. 
Symbols represent the binned luminosity function in bins of ~0. 1 
dex. 

Alongside each panel of LFs, we include a graphic representa- 
tion of the error ellipses of the 1, 2 and 3a contours in L*-(j>* space 
representing the boundaries of the volume in 2D (<jf>*, L*) space. 
We anchor a to the median value across the redshift slices for each 
line as shown in each Figure, while L* and <f>* are allowed to vary 
(i.e. a is constant per panel to allow comparison of LFs). 



on the size of the aperture used. As our analysis is carried out in 2 
arcsecond apertures and a retrospective aperture correction applied 
(see below), we assess completeness in the same manner down to 
the 5a limiting magnitude. We construct a series of images insert- 
ing randomly positioned fake point sources in each image, allow- 
ing us to determine recovery fraction as a function of magnitude. 
We estimate incompleteness per AniNB = 0.05 bin, and create an 
appropriate number of sources with NB excesses distributed fol- 
lowing the observed distribution of detected objects, and with NB 
magnitudes uniformly distributed within the bin. The effect is small 
as we only populate the faintest bins, typically adding fewer than 
10 objects to a LF. This number rises to 21 objects (a 2 per cent in- 
crease) for the second most densely populated LF ([O II] in NB816) 
and 91 objects (a 4 per cent increase) for [O II] in NB921 where the 
number of detections is more than a factor of 2 larger again. Cor- 
rections are particularly low due to the stringent 5a NB -detection 
requirement. 



5.4 Galactic Extinction and Aperture Correction 

Each LF is corr ected for Galactic ext inction of Ay = 0.070 magni- 
tudes using the Card elli et all i 19891) reddening law (following the 
same method as in Section l4~TV As the analysis is carried out using 
aperture photometry, we must assess the fraction of flux that is lost 
by only considering the light within our 2 arcsecond aperture. We 
calculate the ratio of total NB flux to aperture NB flux using the 
output SEXTRACTOR values (FLUX_BEST and FLUX_APER respec- 
tively) for each redshift slice, and scale each output LF according 
to the median ratio in that slice. The correction is typically small 
(a factor of ~ 1.2) even at the lowest redshifts, however for the 
smoothed VISTA NB data we find a slightly larger correction (~ 
1.8) due to a much broader PSF. 



6. 1. 1 Treatment of Shallow Data 

As discussed in Section [5] the maximum likelihood analysis min- 
imizes Equation [8] allowing (/>* , L* and a to vary. However, this 
analysis requires a deep and plentiful data set in order to produce 
a well constrained fit, particularly for a. NB570 is the shallow- 
est of the Subaru NB images, and due to its blue wavelength sam- 
ples [O III] emitters at a very low redshift (in the smallest cosmic 
volume of the analysis); consequently the LFs produced from the 
NB570 filter selection ([O III] at z=0.14 and [O II] at 0.53) are 
subject to low number statistics, with only 42 and 83 objects re- 
spectively making up each LF. 

Furthermore, VISTA data are ~1 mag shallower than NB570, 
and so our data struggles to probe L* in the z = 0.5, 1.0 and 1.6 
redshift slices, in addition to very low numbers of emitters being 
detected. 

The shallower depth of this data means that a is not as well 
constrained as for the deeper more detection-rich LFs, however, 
the NB985 combined data set is selected close in wavelength to 
the NB921 filter, and consequently probes redshift slices only 
marginally more distant than the NB921 selection. To combat the 
shallow data we fit LFs for the 3 VISTA redshift slices with a con- 
strained by a Gaussian prior based on the value for the NB921 se- 
lection for the same emission line. This is a valid method of evalu- 
ation where we do indeed have some prior knowledge of the value 
of a. We take a similar approach with the NB570 [O II] LF at 
2 = 0.53, setting the prior according to the [O II] LF at z — 0.35. 

6. 1.2 The Effect of Cosmic Variance 

In order to study the effect of our limited number of detections we 
undertake a jackknife analysis in every redshift slice. For each slice 
we use random sampling with replacement to construct 1000 syn- 
thetic catalogues of narrow-band emitters such that each sample is 
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of the same size as our true catalogues. We fit luminosity functions 
to each of these synthetic catalogues and measure the standard devi- 
ation in each parameter. These la deviations are comparable to but 
no bigger than the uncertainties in the true dataset. Errors quoted in 
Table|4]refer only to the la uncertainty in the real data. In Figure 
|9]we plot two error bars for each point. The smaller of the two er- 
ror bars represents the Poissonian error, and the larger includes the 
standard deviation on the number of objects in each luminosity bin 
from the jackknifed samples added in quadrature. 

6.2 Star Formation Rates 

We use the best-fitting Schechter function in each of the 12 red- 
shift slices to calculate the characteristic SFR and psfr across 10 
G yr of cosmi c time. SFRs are calculated according to the relations 
of lKennicuttl Il998h for Ha (Equation [Tit. [O II] (Equation [T6t 
and using the standard line r atios: Hq/H/?=2.78 and [OIII]/H/J=3 
dOsterbrock & Ferlandl 120061) to derive Equation [TBI for [O III]. 
E quation 1 151 is in good agreement with the empirical relationship 
of iLv et alTl2007b . This method determines extinction-corrected 
SFRs assuming continuous star formation, Case B recombination 
and a Salpeter IMF. 



results in 7 = 3.02, similar to the parametrization of Sobr al et"al1 
d2013h . The discrepancy between these results and shallower slopes 
in the literature is attributed to the small volumes probed in our low 
redshi ft slices. For instance, using the figures of ISomerville et al] 
(2004) we expect uncertainty in psfr in the lowest redshift slices 
to be comparable to the measurements themselves due to cosmic 
variance. Larger areas of sky are required for this technique to be 
of use in determining the psfr at low redshift as the redshift slices 
sampled in narrow-band surveys are especially susceptible to cos- 
mic variance. We re turn to this po i nt in Section 17.41 We note also 
that the [O II] LF of Gilbank et al.1 d20ld) at z ~ 1 gives a psfr in 
good agreement with the progression of [O II] values seen in panel 
[Mb). 

This technique is however sensitive to particularly low SFRs, 
probing well below SFR* (<s£. IMq yr _1 ) in the deepest data sets 
allowing a more accurate evaluation of a. While the role of these 
faint galaxies' contributions to the overall psfr is unclear, if we 
wish to understand the nature of this faint population it is important 
to first constrain a. 



SFRtMeyi-- 1 ) = 7.9 x 10" 35 L Ha E(Ha) 



(14) 



SF^Meyr" 1 ) = 7.35 x 10" 35 L [OIII] E[OIII] (15) 



SFR(M Q yr- 1 ) = 1.39 x 10" 34 L [O n]E(Ha) 



(16) 



Figure [l0](a) presents the fully-corrected characteristic star 
formation rate (SFR*; aperture and extinction corrections described 
in Sections 16.11 and |5.4t for each redshift slice in this analysis, 
red points showing Ha derived estimates, green, [O III], and blue, 
[O II]. Error bars represent SFRs for the corresponding la upper 
and lower bounds on the characteristic SFR. Figure [Tol b) presents 
the integrated psfr for each LF (fu ll integration , extin ction- 
corrected) together with the fits of ISobral et al.l ( 2013b and 



Hopk ins & Beaconj (2006) using the data o f Baldry & Glazebrook 
d2003l) . 



7 DISCUSSION 

7.1 Evolution of the Star Formation Rate 

We find characteristic SFRs broadly consistent between all three in- 
dicators out to z ~ 1.6. The characteristic SFR increases by a fac- 
tor ~ 10 between z ~ 0.14 and z ~ 1.6, in good agreement with 
literature results. The Ha and [O II] derived star formation rates 
show a distinct increase with redshift consistent with other studies, 
but we note that the [O III] estimates (while broadly consistent with 
the other indicators due to their uncertainty) would produce a more 
convincing increase with redshift if for instance the [O III]/H,9 ratio 
were to decrease with redshift. 

The parametrization of psfr. according to psfr oc (1+z) 7 
gives 7 = 4.58 when handling errors as the maximum and mini- 
mum values of psfr produced from the la confidence limits of 
each measurement of <f>* , L* a nd a. We note t hat applying the same 
least squares fit to the data of ILv et al.l ( 120071) gives a slope of 7 = 
4.20. Applying a blanket 10 per cent error on psfr to our data 



7.2 The Faint End Slope 

The gradient of the faint end slope, a, i s the subject of much debate. 
The recent work of ISobral et al.l J2013) made a significant advance 
on previous studies by examining estimates of a across a broad red- 
shift range using the first sample of consistently selected emission 
line objects out to z — 2.23. An important part of determining the 
most likely luminosity function at any redshift however, is to de- 
termine the significance of the galaxies we are capable of detecting 
in the analysis within the general population at redshift z. As dis- 
cussed in Section [5T2l we compute the ratio of detectable galaxies 
to the total number of galaxies capable of hosting the line. Inves- 
tigating this ratio we find that the best fit a of each LF is in fact 
very sensitive to the assumed faintest broad-band magnitude of a 
galaxy that can host an em ission line (i . e. th e denominator of this 
fraction). The analysis of ISobral et al.l d2013l) assesses emission- 
line detection fraction down to the limit of their NB data, and so 
we investigate the effect of assessing our data in the same way. 

Evaluating our data set with mfaint set to the limiting NB 
magnitude (mf a j nt (NB)) we find that the resultant values of a are 
consistently steeper than for our mf a i nt (EW) analysis, the medi- 
ans of the two samples being offset by 0.3. Furthermore, the range 
in values of a is reduced (from 0.90 to 0.66) through use of the 
(wifaint(NB)) limit across the same redshift range. 

Taking the well-constrained NB921 Ha and NB921 [O II] 
luminosity functions we vary the input limiting equivalent width 
value for mf a i nt . The sensitivity of a to the input value is displayed 
in Figure QT| As EW mfaint is increased, a larger fraction of the 
population of galaxies can host emission lines, and so the detection 
fraction decreases and the resultant a steepens. Conversely how- 
ever, as the input EW mfaint is lowered below 100 A the numerator 
of Equation [T2] (the red portion of a line in Figure [8} begins to de- 
crease when mfaint crosses the 5a detection limit. As the detection 
limit moves to brighter magnitudes, the galaxies undetected due 
to their position beneath the horizontal observed equivalent width 
limit in Figure [8]become a larger fraction of the overall population, 
and so the detection fraction again decreases, and a steepens. The 
shallowest a therefore, should correspond to the point where the 
NB detection limit, ms„, is equal to mfaint. 
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Figure 9. Graphic representation of our best-fitting luminosity functions (Gaussian prior on a for VISTA LFs and NB570 at z = 0.53) and their 1, 2 and 3<r 
error ellipses (a fixed to the median value per line for comparison) in L*-<j>* space. The smaller of the two error bars on each point in the lefthand panels are 
Poissonian, larger error bars have the la variation as estimated by our jackknifing technique added in quadrature. Ha (top panels) [O III] (central panels) and 
[O II] (lower panels). The upper SFR scale on these plots gives an indication of the raw SFR values for these LFs, these numbers have not been corrected for 
dust extinction. 
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Table 4. Luminosity functions in 12 redshift slices for Ha, [O III] or [O II] emission. Numbers of objects quoted refer to the final number of objects per 
luminosity function (after culling to the broadest 3<r colour-scatter limit of the SXDS fields, and taking into account the completeness correction). Volumes 
refer to the volume enclosed where the filter transmission is greater than 50 per cent. Schechter fits to the data are presented with the value of L* corrected for 
aperture effects and Galactic extinction at the wavelength of the filter. The extinction corrected values listed are those derived for 1 mag of extinction at Ha. 
Where the raw data did not produce a converging LF, values are absent for SFR* and Psfr as these would be non-physical. 
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Figure 10. Panel (a) shows the extinction-corrected characteristic star formation rate (SFR *) determined in each redshift slice, error bars corresponding to the 
f<r error on L*. Colours reference the emission line used to derive the SFR. Panel (b) presents p sfr in each redshift slice. Colours again represent emission 
line used to derive the measurement. The solid black line shows our parametrization of p sfr oc (1 + z) 4,58 , and the dotted lines show parametrisations in 
the literature. 



7.3 The Effect of a on psfr 

The effect of a on psfr is dependent on the value of a itself. For 
shallow a, the effect is small, however, for steeper slopes the psfr 
estimate is altered more dramatically for a small change in a. For 
instance, taking the well-constrained NB921 [O II] LF at z = 1.46 



and altering a = —0.91 by Aa = 0.1 in either direction (hold- 
ing 4>* , L* constant) causes the psfr to rise by ~ 5 per cent for 
steeper (more negative) a but only decrease by ~ 3 per cent for 
the corresponding Aa = 0.1 towards shallower (less negative) a. 
Taking the same approach with the NB921 Ha LF at 2 = 0.4, 
a — — 1.14 , and the equivalent Aa = 0.1 changes in a result in a 
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Figure 11. The variation of a and Psfr with differing input values of the 
limiting equivalent width for mj a ; nt . The upper panel shows the variation in 
resultant values of a for the NB921 [O II] luminosity function at z = 1.46 
(blue) and the NB921 Ha luminosity function at z = 0.40 (red) when the 
limiting value of equivalent width for mf a ; nt is varied. The lower panel 
shows the corresponding resultant values of Psfr as m faint changes for 
the same two luminosity functions. 



i—K colour (psfr = 0.0138 M©yr 1 ). The recovery of this psfr 
despi te the differing aperture sizes used for detection (Sobral et al. 
120131 use a 3 arcsecond aperture at low redshift) lends support to 
the success of our global corrections for aperture flux loss. 

As discussed in Section [6.1.2l the uncertainty in our results due 
to the finite number of objects detected per redshift slice is inves- 
tigated via a jackknife analysis. However, it is not possible to tell 
if the values of psfr we compute are representative of the overall 
Psfr at each redshift. This depends on the large scale structure in 
the UDS at the redshifts we sample. It is not possible to tell whether 
these regions are particularly under- or over-dense without a com- 
parison to many more similar surveys. For example, in the lowest 
redshift slice that we sample, [O III] at z — 0.14, we would require 
a contiguous volume ~ 140 times larger in order for the la uncer- 
tainty in number counts of galaxies to reduce to 2 per cent accord- 
ing to the calculations of lSomerville et al. (2004). The issue can be 
more efficiently addressed by observing widely separated regions 
of sky and comparing the variations seen to between such fields to 
predictions. For exam ple the divergence o f psfr by a factor of ~ 3 
between this work and lSobral eUtud2012l) 's UDS+COSMOS mea- 
surement at low redshift (Ha at z = 0.40) is quantitativ ely con- 
sistent with the predictions of of ISomerville et all d2004l) . Finally 
we note that our value of p sfr at z = 1 .46 is in excellent agree- 
ment with both th e p SF R oijSobral et all d2012l) 's UDS+COSMOS 
measurement, and lLv et"aL I d2007h 's Subaru Deep Field (SDF) mea- 
surement. Based on these comparisons we note that our estimates 
of Psfr in the lower redshift slices may differ from the true values 
by up to a factor of ~ 8 in the lowest redshift slice, z — 0.14. 
The resultant steep value of 7 should therefore not be a cause for 
concern when the uncertainty due to cosmic variance is so large. 



1 1 per cent increase in psfr for steeper a and a 8 per cent decrease 
for a shallower slope. This behaviour is reflected in the lower panel 
of Figure[TT] however it is also apparent that psfr increases more 
rapidly towards smaller limiting EW values for mfaint- We inter- 
pret this as the result of an unrealistic limit to mfaint causing the 
detection fraction of galaxies to become very small and produce an 
unrealistic correction inconsistent with the data. 

7.4 Comparison to Other Studies 

The NB921-selected [O II] LF and p SF R at 2 = 1.46 is 
in excellent agreement with the results of ISobral et al l d2012l) 's 
UDS+Cosmic Evolution Survey (COSMOS) selection. Conversely 
though, comparison of our NB921-selected sample of Ha emit- 
ters to the selection of Sobral et al. (2013) at z = 0.4 reveals 
our psfr is much l ower than the UDS+COSMOS value. However, 
Sobral et all d2013l) found a difference in psfr between their UDS 
and COSMOS-selected samples, with the UDS value being a factor 
of 3 lower than that for COSMOS. Their best-fitting LF using UDS 
emitters alone produces psfr = 0.014 Moyr - (D. Sobral, pri- 
vate communication) a value slightly higher than our own. We are 
able to fully account for this diffe rence if we use the s ame method 
of redshift slice assignment as in ISobral et al. I d2013l) (where ob- 
jects with K mag > 23.0 were assigned according to B—R, i—K 
colour). Our work uses a later UDS data release and so deeper near- 
infrared imaging results in more accurate photom etric redshifts for 
faint o bjects — these are in good agreement with Griitzbauc h et alj 
d201ll) . We also confirm the reliability of our apertur e corrections 
throug h the agreement of psfr at z = 0.4 between 1 Sobral et"al] 
d2013l) 's UDS-only sample and our Ha emitters according to B—R, 



8 CONCLUSIONS 

We have developed a robust method for evaluating luminosity func- 
tions which correctly handles the selection criteria of NB samples, 
negating the need for retrospective corrections to the data set. As 
a result, the method is suitable for highly irregular filters, and is 
applicable to a broad range of data sets. This is the first maximum 
likelihood analysis of narrow-band emitters using one of the largest 
NB-selected data sets to date. Our results can be broadly summa- 
rized as follows: 

• We have presented maximum-likelihood luminosity functions 
in 12 redshift slices out to z = 1.64, and demonstrated that this 
technique produces results in concordance with literature luminos- 
ity functions using similar data sets ove r similar redshift ra nges. 
In particular we c onfirm the results of ISobral et al.l d2013l) and 
ISobral et all d2012l) in the z = 0.4 and z = 1.46 redshift slices 
for the Ha and [O II] NB921 samples respectively. 

• We find the evolution of the psfr to rise according to psfr 
oc (1+z) 4 ' 58 but conclude that this is strongly affected by cosmic 
variance in the small volumes probed at low redshift through this 
technique. 

• While the areal coverage of this survey is not large enough 
to place constraints on psfr at low redshift, the sensitivity to very 
low SFRs that deep narrow-band observations provide makes this 
technique ideal for constraining the faint end slope. 

• We discover that the detection fraction of emission lines plays 
an important role in the determination of the faint end slope. We 
propose that the limit of integration for the total number of galax- 
ies capable of hosting an emission line should be empirically mo- 
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tivated, and here we base this on the line of interest's maximum 
observed equivalent width in the local Universe. 

• The deep NB data will still provide useful information on the 
population of objects making up the faint end slope, and allow an 
evaluation of SFR with stellar mass. We will present this analysis 
in a future paper. 
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